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I.  INTRODUCTION 

In  discussions  of  chaff  cloud  tracking  it  is  often  tacitly  assumed 
that  the  measured  radar-cross-section  (RCS)  centroid  follows  a  Keplerian 
orbit.  This  is  not  entirely  true  for  two  reasons: 

1.  RCS  fluctuations  cause  the  instantaneous  cross-section  centroid 
to  move  randomly  about  the  long-term-average  centroid.  This 
means  that  between  any  two  independent  measurements  there  will 
be  random  fluctuations  in  the  centroid  position,  even  with  a 
noiseless  radar . 

2.  For  a  large  cloud  there  are  various  nonlinear  gravity  gradient 
effects.  For  typical  trajectories  the  long-term-average  cen¬ 
troid  of  a  cloud  will  drift  from  a  Keplerian  particle  having 
the  same  initial  velocity  and  position.  This  displacement 
distance  is  quite  small,  however,  being  on  the  order  of  0.1% 
of  the  cloud  radius  even  for  rather  large  clouds. 

Each  of  these  effects  is  discussed  in  subsequent  subsections. 
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II.  CENTROID  MOTION 

A.  ANALYSIS 

Consider  a  cloud  that  is  N  range  resolution  cells  long  but  small 
in  angular  extent.  The  centroid  is  defined  as 

N  /  N 

C  =  ]Ci0i  /£  °i  h  y/z  (1) 

i=l  /  i=l 

where  o  is  the  RCS  of  the  ith  cell. 

If  all  of  the  c(i)  were  equal  and  constant  in  time,  then  we  would 

have 

=  N(N+l)/2  N+l 
C  N  "  2 

Thus  if  there  were  25  resolution  cells  the  centroid  would  be  in  the 
center  one;  that  is,  (25+l)/2  =  13  . 


Because  of  the  nature  of  chaff  responses,  however,  the  RCS  will  not 
be  constant  from  cell  to  cell  but  will  fluctuate  because  of  interference 
effects.  If  the  number  cf  chaff  pieces  in  any  cell  is  sufficiently  large 
the  probability  density  function  p^(o)  for  the  distribution  of  indepen¬ 
dent  measurements  of  one  cell  will  be  exponential;  that  is 


Pj_(o) 


1  -o/o 
—  e 


o 


We  shall  for  the  time  being  consider  the  case  in  which  the  expected  value 
of  a  is  the  same  for  each  cell. 
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The  centroid  of  the  first  N  cells  can  be  written 


c(N)  = 


N  N-l 

E  ioi  Non  +  E  iai 

i=i  _  w  j=i  1 

N  =  N 

E  o±  E  ^ 

t  _  1  A  *  1  i 


which  can  be  rearranged  to  give 


(n)  E° i 


NaM  +  c (N- 


-»  £  °i 


If  we  let  gN  =  c(N) 


(#,  ••) 


then  we  have  the  contiguity  relationship 


gN  %-l  +  N°N 


The  probability  of  having  a  value  g^  is 

'V  -  /  f'sn-i)''! (“ii  -  —  r-1)  d%-l 


Taking  Laplace  transforms  we  have 


V>  =  W> 


1  +  Nsa 


where  we  have  used  the  face  that 


i  \  1  -o/a 

Px(a)  =  -3-  e 

a 
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Similarly 


PN-1(S)  "  PN-2(S) 


1 

1  +  (N-l)so 


and  thus 


TV-1 


and 


But 


PN(s)  -  Tf  u  +  iso) 

i=l 


N 

£nP.,(s)  -  -  j*  £n(l  +  iso) 
i-1 


N 


^  ^  £n(l  +  iso) 
i-1 


- .  iHii  +  Mi 

2  3 


.  r.7-i 

[• 


- 1 


+  |(2s7)  .  +  iz&l . .  .  . 


+  j(Ns7)  .  iSKli  +  jfflsati  .  .  .  .] 

—  2 

-  (so)  [1  +  2  +  3  •  •  •  N]  -  [l2  +  22  +  32  •  •  •  N2] 

+  ('3^>-  [l3  +  23  +  33  •  •  •  N3]  -  .  .  . 

=  As  +  Bs2  +  Cs3  .  .  .  =  R(s) 
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.at  n-  - 


where  A  =  aN(N+l)/? 

B  =  “I”  [N(N+1)  ■^"^'' 

a3  9 

C  =  j-  [N(N+1)/2]Z 


P<%>  =/"\^)  V^"1  ^"H(S) 


Since  g  =  (c  )  £  a 

i-1  1 


we  can  write 


P(g)  =  /  P(c)PN(z  =  g/c)  ~  dc 


where  we  have  suppressed  the  subscript  on  g  and  c  ,  and  let 


a.  ,  z 


Before  going  on,  it  is  worthwhile  pointing  out  explicitly  that  in 
Eq.  2  we  have  tacitly  assumed  that  p(c)  was  independent  of  z  .  Consi¬ 
der  two  sets  of  o's,  a^,  o2,...a^  and  a 0^,  ,  each  having 

the  same  total  o  (i.e.,  °i+a2’  * '°N  =  °l+a2’**°N  =  z  )•  Each  set  has  a 

unique  centroid  (c  and  c)  associated  with  it.  The  relative  likeliness 
of  each  set  (point  in  N  dimensional  space)  is 

p(o;1)p(a2) .  ..p(oN)  da^  dc?2...daN 

p(o1)p(o2) . .  .p(aN,  do1da2...doN 


Now  consider  two  other  sets  of  points,  ota^.aOj, . .  .aa^  and 
aa^jaa^, • • .aa^  .  These  have  the  same  centroids  (c  and  c)  associated 
with  them.  If  the  relative  likeliness  of  these  two  new  sets  of  points 
is  the  same  as  for  the  first  two  sets  (a  *  1),  then  the  density  dis¬ 
tribution  of  c's  will  be  invariant  with  a  .  Since  varying  a  varies 
the  total  RCS  (z) ,  we  will  have  proved  that  p(c)  is  independent  of 
z  . 


Thus  we  must  show  that 

pCo^  . .  .p(oN)  do^.  •  -daN  ,  p(aa1>  . .  .p(aoN)  dac^. .  .daaN 
p^)  . .  .p(aN)  do^ . . . do^  p  (aa^)  . .  .p (aa^)  d(aa1)  . .  .d(aaN) 

Canceling  the  da's  and  the  a's  gives 

pCo^^)  . .  .p(aN)  ?  p(aa^) . .  .p(aa^) 

P^)  . . -p(oN)  p(aa1) . .  .p(aaN) 


For  an  exponential  distribution  this  gives 


-z  „  -az 
e  ?  e 


-z  -az 
e  e 


and  the  equality  has  been  shown. 


The  expected  distribution  of  J]  o,  ,  i.e.,  pN(z)  »  can  be 

i=l  1  W 

evaluated  by  noting  that 


£ 

3N(Z)  =  /  pN- 


1(o')p1(z  -  o')  do' 


The  N  indicates  a  summation  over  N  cells.  Taking  Laplace  transforms 
we  obtain 


r 


b> 


i 


i 


It  K 


■ 


i 


KJ 


PN(S)  =  tPN-l(s)]Pl(s) 


[PN-2(S)P1(S)1P1(S) 


[P^s)]1 


But 


P^s)  =vf~  e"0/o 


a 

1 


1  +  as 


Thus 


V*>  vf 


-1 


(1  +  as)^ 


.  IsZs)?  1  e-*/° 

a"r  (N) 


If  we  multiply  both  sides  of  Eq.  2  by  gq  and  integrate  from  zero 
to  infinity,  we  have 


CO 

/ 


cq  p(c)  dq  = 


TN 


71  r (N+q)  ^ 


00 

/ 


gq  p(g)  dg 


If  we  let  q  =  1  ,  then  we  have 


The  integral  on  the  left  is  c  (the  expected  value  of  the  centroid). 
The  integral  on  the  right  is  the  first  moment  of  p(g)  and  this  is 
simply 


e-H(s) 


s=0 


=  A 


Thus 


,  /H-  N+l 

C  =  A/N 0  =  — r— 


Similarly,  if  we  take  q  =  2  ,  we  have 


00 

/ 


c  p(c)  dc 


00 

f 

o2(N+1)!  J 


g  P(g)  dg 


Remember  that  if  F(s )  =  /  e  SXf(x)  dx  ,  then  differentiating  both 

'0 

sides  N  times  with  respect  to  s  gives 


(if  F<s>  -  <-1)S  /  £« 


N  -sx  , 
x  e  dx 


If  we  now  set  s  =  0  ,  we  have 


1 xNf(x)  dx  =  (_1)N  ih)  F(s) 


s=0 


>N 


Thus  the  Nth  moment  of  a  function  is  (-1)  times  the  Nth  derivation  of 
its  transform  at  s  =  0  . 


"2?  J&&**t  t 


0  (N+1)N 


(A  -  2B) 


The  variance  is 


—  2  2  2 
2  -  A  -  2B  A 

c  ~  c  =  — 2  '-22 

a  N(N+1)  oN 


In  order  to  verify  this  result  a  computer  Monte  Carlo  experiment 

was  run.  The  quantity 
N 

2>i 
i=l  x 
c  =  __ - 

E  *± 

i=l 

was  calculated  after  generating  N  random  values  of  .  The  o^  were 
selected  from  a  distribution  of  the  form 


P(a)  =  e  ° 

The  process  was  then  repeated  W  times  to  give  a  total  of  W 
independent  values  of  c  . 

A  total  of  W  =  500  independent  evaluations  of  c  ,  for  a  cloud 

containing  121  cells,  gave  a  value  of  c  =  61.1  ,  versus  a  theoretical 

2  2 

value  of  61,  i.e.,  (121  +  l)/2  .  The  quantity  c  -  f  had  a  value  of 

121  -  1 

10.1  compared  with  the  theoretical  value  of  10.0  (i.e.,  - — —  )  • 


w? .3  vs.v  *■«  ~ 
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B.  INTERPRETATION 

Since  the  variance  is  proportional  to  N  (for  large  N  ) ,  the  rms 

jlf 

fluctuation  (standard  deviation)  on  any  measurement  is  proportional  to 
v^j  .  If  the  cloud  has  a  total  length  L  and  a  resolution  cell  width 
6  =  L/N  ,  we  see  that  on  any  measurement  the  instantaneous  "true"  centroid 
is  typically 


l  =  standard  deviation  *  ceil  size 


(3) 


from  the  long-term-average  centroid. 


This  result  is  physically  reasonable  in  that  the  fluctuation  gets 
smaller  as 


1.  The  cloud  gets  shorter 

kk 

2.  The  resolution  gets  better 


The  "square  root  of  N"  form  for  the  deviation  from  center  is  exactly 
what  one  would  expect  from  a  random  walk  configuration  and  could  probably 
have  been  predicted  on  a  priori  grounds;  however,  the  coefficient  of  1/12 
could  not  have  been  similarly  predicted. 


It  would  be  slightly  misleading  to  refer  to  this  term  as  "error"  since 
it  is  the  true  instantaneous  RCS  centroid.  The  point  is  that  the  RCS 
centroid  fluctuates.  The  most  one  can  say  is  that  there  is  some  long¬ 
term-average  centroid  which  can  be  used  as  a  reference  point.  This 
point  is  not  necessarily  coincident  with  the  center  of  mass.  Strictly 
speaking  one  can  argue  that  even  this  long-term-average  has  a  secular 
variation  because  of  changing  asp'  t  angle  of  the  radar  line-of -sight 
relative  to  the  chaff  spin  axes.  This  effect  is  small. 

Remember  we  have  assumed  there  are  a  number  of  chaff  pieces  in  each  cell 
(to  assume  an  exponential  distribution  of  RCS  values)  and  also  that 
L>>o  (N>>1)  .  These  two  assumptions  impose  limits  on  where  the  result 
can  be  used. 
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This  brings  up  the  interesting  point  as  to  whether  there  might  be 

an  alternate  partitioning  scheme  which  would  have  smaller  fluctuations 

associated  with  it.  We  have  been  looking  for  the  centroid  which  is  a 
* 

quantity  defined  by 


dx  (x  -  c)  a(x)  dx  =  0 


(A) 


It  is  also  possible  to  track  relative  to  the  median;  one  looks  for 
the  point  that  has  equal  amounts  of  RCS  on  either  side. 


Mathematically  we  look  for  a  point  m  such  that 
m 


J  a(x)  dx  =  J 

_oo  in 


a (x)  dx 


or,  put  into  the  same  format  as  Eq.  4 


CO 

/ 


x  -  m 

I  x  -ml 


a(x)  dx  =  0 


(5) 


(6) 


The  difference  between  the  two  quantities  is  that  the  calculation 
of  c  places  more  weight  on  pieces  of  chaff  far  removed  from  the  center 
(the  weighting  factor  has  a  magnitude  |x  -  c|  ).  The  definition  of 
m  places  equal  weight  on  all  pieces  of  chaff  since  the  weighting  factor 
in  Eq.  6  always  has  unit  magnitude. 

The  Appendix  gives  a  derivation  of  the  rms  fluctuation  when  the 
median  is  tracked  and  it  is 

l  -  x  6  =  -|  v*LS 


For  generality  we  shall  use  an  integration  notation;  it  is,  mutatis 
mutandis,  consistent  with  Eq.  1. 
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Note  that  the  form  is  the  same  as  for  centroid  tracking 

although  the  coefficient  is  slightly  larger. 

C.  GENERALIZATIONS  TO  OTHER  GEOMETRIES 


1.  Non-Uniform  Density 

Throughout  this  entire  discussion  we  have  assumed  that  the  expected 
chaff  distribution  is  uniform.  If  it  is  not  uniform  but  (as  is  more 
likely)  peaked  in  the  center,  the  variance  is  be  less. 


A  chaff  density  distribution  as  is  shown  in  Fig.  1  would  have  a 
variance  of  approximately  /l/12  /L/2  rather  than  /l/12  /l  .  (The 
wings  are  assumed  to  be  negligible  compared  to  the  center.) 


12 


•  *  3P  *,-Y 
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If  the  density  distribution  were  as  shown  in  Fig.  2,  the  centroid 
variance  could  be  considered  as  from  two  contributors,  one  of  RCS  p 
and  variance  G(L2>  and  the  other  RCS  p^  and  variance  G(L  )  . 

G(L)  is  the  Green  function  for  the  variance  of  a  uniform  density 
distribution  of  length  L  and  resolution  length  6  .  Thus  G(L)  =  (1/12) 6L 
for  centroid  mapping  and  (3/12) 6L  for  median  mapping. 


The  resultant  net  variance  for  the  two  is 

a2  -  p°L2G^L2^  +  poLlG^  / 
poL2  +  PJ~1 


^•^8^^®  2*  Chaff  Density  Versus  Range 
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where  g  =  1/12  for  centroid  mapping  and  3/12  for  median  mapping.  This 
result  is  the  same  as  that  obtained  before. 

A  triangular  distribution  would  give 


Similarly  a  quadratic  distribution  would  give 

’  -  (l  -  (mf) 

2  3  T. 

a  =  j  gL6 

All  of  the  preceding  results  can  be  summarized  by  saying  that  the 
rms  shift  of  any  measurement  from  its  long  term  average  is 

(✓gf  )  JU 

where  g  =  1/12  for  centroid  mapping  and  3/12  for  median  mapping;  f 
is  a  form  factor  which  is  1  for  a  uniform  chaff  distribution,  3/4 
for  a  quadratic  distribution  and  2/3  for  a  triangular  one.  Thus  for 
example,  the  centroid  of  a  chaff  cloud  of  length  L  with  a  quadratic 
distribution  of  chaff  will  fluctuate  about  its  long-term  position  with 
a  typical  deviation  of  / (1/12)  (3/4)  *^L6"  =  0.25  y^L6  . 


2. 


Three-Dimensional  Configurations 
The  discussion  throughout  this  paper  has  been  restricted  solely  to 
finding  the  variance  of  the  RCS  centroid  (or  median)  of  a  cloud  that  was 
long  in  one  dimension  but  short  in  the  other  two.  (We  have  considered 
the  long  dimension  to  be  range,  but  it  could  just  as  well  have  been  one 
of  the  angular  dimensions.) 


It  is  worthwhile  generalizing  our  result  to  the  case  in  which  the 
object  could  have  subdivisions  in  angle  as  well  as  range. 


Consider  a  rectangular  cloud  (Fig.  4)  of  dimensions  L^,  and  L^, 
which  can  be  radar  subdivided  into  cells  of  size  and  A^  (L^ 

could  be  the  range  extent  of  the  cloud  and  A^  the  range  resolution;  L ^ 
could  be  the  width  of  the  cloud  in  one  dimension,  and  A^  the  associated 
angular  resolution) . 


Figure  4.  Graphic  Representation  of  Chaff  Cloud 
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There  are  thus 


columns  of  length  L^  ,  and  each  one  has  a  centroid  standard  deviation 
of 

JE 

*12 

If  we  average  the  centroids  of  each  of  the  N  columns,  the  stan¬ 
dard  deviation  will  be  reduced  by  a  factor  1/v/’n^~  ,  to  give  a  value 


The  standard  deviation  in  the  other  dimensions  can  be  written 
similarly,  i.e., 


The  standard  deviation  of  the  net  (three-dimensional)  "error"  is 
the  square  root  of  the  sum  of  the  square  of  the  standard  deviation  in 
each  dimension.  Thus  we  have  a  net  standard  deviation  of 


If  we  re-express  the  cloud  in  terms  of  two  form  factors  =  L^/ 
and  a  =  L»/L1  ,  and  one  scale  size  L,  ,  then  we  have 


standard  deviation  =  -  /l/12  %Ja7 aTaT 

fir  >  1  2  3 


9  2 

!+a2  +  a3 


Note  that  for  a  given  shape  cloud  ( a 2  and  constant)  the  standard 
deviation  gets  smaller  as  the  cloud  gets  larger.  This  is  contradistinction 
to  the  case  of  a  one-dimensional  cloud  where  the  standard  deviation  gets 
larger  as  the  cloud  gets  longer.  The  reason  is  that  the  variance  in  any 
one  dimension  is  proportional  to  the  length  in  that  dimension  divided  by 
the  number  of  columns  in  that  dimension,  e.g., 


•  (42a3> 


If  all  dimensions  are  doubled  then  the  variance  of  one  column  is 
doubled  (because  it  is  twice  as  long) :  but  there  are  now  four  times 
as  many  columns  for  a  net  improvement  of  a  factor  of  two  (a  factor  of 
the  square  root  of  two  in  standard  deviation) . 

This  analysis  has  several  limits  which  should  be  spelled  out. 

In  the  first  place,  we  have  tacitly  assumed  no  errors  introduced  by 
S/N  considerations,  i.e.,  the  signal  is  so  strong  that  this  error  is 
negligible  compared  to  variations  in  centroid  caused  by  the  RCS  fluc¬ 
tuations  . 


As  the  cloud  gets  larger  this  assumption  will  become  poorer  for 
two  reasons : 

1.  For  a  constant  amount  of  chaff  (constant  total  average  RCS) 
the  signal  in  each  cell  becomes  smaller  as  the  chaff  is 
diluted,  and  thus  the  signal  noise  error  increases 
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2.  The  centroid  fluctuation  "error"  decreases  (yl/AL). 
Thus  at  some  size  S/N  considerations  will  dominate. 


There  is  a  further  limit.  If  we  consider  the  chaff  cloud  to  be 
compressed  to  a  point  (still  maintaining  a  constant  total  RCS)  then  very 
accurate  range  measurements  are  possible  because  of  the  large  S/N  ratio. 
The  error  will  be  proportional  to 

resolution  length 

/s/N 

and  will  not  have  a  centroid  RCS  fluctuation  term. 


The  competing  effects  can  be  illustrated  as  shown  in  Fig.  5,  where 

we  have  considered  a  cubical  cloud  of  side  L  .  For  small  L  (much  less 

1  *  2 

than  a  resolution  length  and  even  less  than  resolution  length/ (S/N) 
the  centroid  standard  deviation  comes  solely  from  S/N  and  radar  consider¬ 
ations  and  not  from  true  motion  of  the  centroid  (Region  I) . 


As  L  gets  larger  (Region  II)  centroid  motion  comes  into  dominance 
and  the  standard  deviation  becomes  larger.  After  L  becomes  larger 
than  a  resolution  cell  and  the  radar  has  multiple  cells,  the  centroid 
motion  error  becomes  smaller  (at  a  rate  proportional  1/A.)  (Region  III). 


Finally  (Region  IV)  the  standard  deviation  caused  by  the  centroid 
motion  becomes  less  than  the  increasing  noise  error  (because  of  decreased 
S/N  in  each  cell).  D.  Hunt  has  shown  that  in  this  region  the  standard 
deviation  of  the  centroid  error  is  independent  of  cloud  size,  as  long 
as  the  S/N  ratio  is  large  enough  to  permit  detection. 

jKj 

D.  Hunt,  Chaff  Cloud  Centroid  Measurements  with  the  MSR  and  ALOR  Radar 
Systems,  General  Research  Corporation  IMR-1377,  August  21,  1970. 
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STANDARD  DEVIATION  OF  CENTROID 


AN-26659 


III.  GRAVITY  GRADIENT  EFFECTS 

In  discussions  of  chaff  tracking  it  is  usually  assumed  that  the 
(long-term-average)  centroid  follows  a  Keplerian  orbit,  although  it  is 
known  that  this  is  not  rigorously  true.  To  quantify  the  magnitude  of 
the  error  in  such  an  assumption,  several  machine  runs  were  performed. 

The  first  was  a  4500  n  mi  reference  trajectory  30  deg  reentry  angle. 

Then  the  trajectories  of  two  particles  ejected  53  s  after  launch  and  in 
opposite  directions  to  each  other,  were  calculated.  The  center  point 
of  the  two  ejected  particles  can  be  considered  as  the  centroid  of  a 
rather  degenerate  (2  piece)  chaff  cloud.  The  deviation  of  this  cen¬ 
troid  from  the  reference  trajectory  is  a  measure  of  its  deviation  from 
the  Kepler  orbit  since  both  trajectories  had  the  same  position  and 
velocity  at  ejection.  (One  point  in  position-velocity  space  is  suffi¬ 
cient  to  define  a  Keplerian  trajectory.) 

Figure  6  shows  the  separation  between  the  centroid  and  the  true 

Keplerian  position  for  four  ejection  velocities  (10,  50,  100,  and  150  ft/s) 

as  a  function  of  time  after  ejection;  the  particles  are  ejected  at 

right  angles  to  the  velocity  vector  and  in  the  plane  of  the  orbit.  Note 

3 

that  the  slope  of  the  curve  is  approximately  proportional  to  t  and 

thus  the  centroid  has  a  velocity  away  from  the  reference  trajectory  at 

2 

a  rate  proportional  to  t  .  The  curves  can  be  very  closely  approximated 
by 

-12  2  3 

D  =  5  x  10  V  tJ  ft 
o 

Figure  7  shows  the  fractional  displacement  of  the  centroid,  that 
is,  the  displacement  divided  by  the  cloud  diameter,  as  a  function  of 
diameter  at  impact.  The  fractional  rms  is  quite  small  and  even  for  a 
cloud  as  large  as  100  n  mi  at  reentry  is  less  than  0.2%.  This  deviation 
is  usually  negligible  for  most  applications. 

* 

We  are  assuming  here  that  each  piece  has  the  same  RCS. 
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10'4  10'3  10'2 
FRACTIONAL  ERROR 

Figure  7.  Fractional  Displacement  of  Centroid  as  Function  of  Cloud 
Diameter  at  Impact:  30-deg  reentry;  4500  n  mi  Trajectory 

If  the  particles  are  ejected  along  different  axes,  the  same  general 
conclusions  hold,  although  numerical  values  are  slightly  different 
(Fig.  8).  The  largest  fractional  shift  of  the  centroid  occurs  when  the 
particles  are  ejected  perpendicular  to  the  plane  of  the  trajectory.  Even 
in  this  case  the  fractional  deviation  is  less  than  1%  even  for  a  ICO  n  mi 
diameter  cloud.  This  is  the  deviation  for  the  special  case  in  which  all 
of  the  particles  are  at  the  edge  of  the  cloud;  any  other  weighting  would 
reduce  the  values  still  further.  A  uniform  distribution  reduces  the 
fractional  error  by  a  factor  of  3.  As  an  example,  chaff  ejected  perpen¬ 
dicular  to  the  plane  of  the  trajectory  (worst  case)  so  that  at  reentry 
it  formed  a  uniform  line  40  n  mi  long  would  have  a  deviation  of  less 
than  0.1%  of  40  n  mi. 

And  finally  a  realistic  case  in  which  particles  were  ejected  in 
several  azimuths  would  give  a  further  cancellation  because  of  offsets 
in  different  directions . 

In  short,  this  gravity  gradient  effect  is — for  typical  trajectories 
negligible. 
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DIAMETER  AT  IMPACT 


Figure  8.  Fractional  Displacement  of  Centroid  as  Function  of  Cloud 
Diameter  at  Impact:  30-deg  Reentry;  4500  n  mi  Trajectory 
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IV.  CONCLUSIONS 

The  major  "error"  in  centroid  tracking  (over  and  above  radar  bias 
and  S/N  errors)  arises  from  the  fact  that  the  centroid  is  not  a  steady 
quantity  but  rather  fluctuates  in  space.  The  rms  fluctuation  from  mea¬ 
surement  to  measurement  is  on  the  order  of  0.3/L6  in  range. 


The  error  caused  by  the  drift  of  the  centroid  (or  median)  from  a 
Kepleiian  orbit  is  small  in  comparison  with  this  effect.  (A  0.1%  frac¬ 
tional  shift  in  the  centroid  would  be  comparable  to  the  centroid  jitter 
effect  only  if  vft  were  on  the  order  of  10  and  this  would  demand  a 
resolution  cell  width  one-millionth  of  the  cloud  length.) 


APPENDIX  A 


In  this  appendix  we  shall  estimate  the  variance  of  the  median  RCS 
position.  The  cloud  is  N  =  L/6  resolution  cells  long,  where  L  is 
the  cloud  length  and  6  is  resolution  cell  width.  Each  cell  has  the 
same  average  RCS  a  ,  and  independent  values  are  distributed  exponen¬ 
tially. 


If  one  sums  the  RCS  values  in  the  first  v  cells  the  probability 
distribution  for  the  sum  is 


(o,/o)V  1  ,  -o./o 

pi(orv) - rw~~' 

0 


(7) 


while  the  probability  distribution  for  the  sum  of  the  remainiv  N-v 
is 

,  vN-(v-l) 

,  ,  (02/0)  1  -°2/a 
P2(°2’V)  =  - - 7  6 


(8) 


If  z  =  , 


then 


p(z 


00 

»v)  =  J  p^a. 


v)  p2(a-z,v)  da 
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For  z  =  0  (equal  RCS  on  each  side) ,  we  have 


00 

>,v)  =  f 


pp(a,v)  p2(a,v)  da 


This  gives  the  distribution  of  the  cell  number  (v)  at  which  the 
RCS  sums  on  each  side  are  equal;  we  shall  name  the  v  value  which 
partitions  the  cloud  into  two  equal  parts  m  and  write 


P (0,v)  =  p(m) 


00 

"  f  IT1 

J  or  (i 


,  s  m-1  1 

(a/a)  - 


.  xN-m-1  -2a /a  da 
(a/a)  e 


of  (m) 


or(N-m) 


1 _ r(N-l) 

a  2N_1  T(m)  r(N-m) 


For  large  N  we  can  make  the  transition  from  a  summation  to  an  integral 
and  write 


A1! 

-I 


dm  m  p(m)  /  /  p(m)  dm 


N/2  /  N/2 

=  J"  dm(nrl-N/2)  p(nri-N/2)  ^  dm  p(m+-N/2) 
-N/2  /  -N/2 


f  dm  (m+N/2)  r(m+N/2)  r (N/2  -  m) 

_  -N/2 _ 

N/2  . . 

f  dm  T (m+N/2)  T(N/2  -  m) 

-W2 

N/2  1 

N  i/2  ^  m  r^/2  +  m)  r(N/2  -  m) 

=  _  +  - 

/*  _ ora _ 

,»  r (N/2  +  m)  r(N/2  -  m) 
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The  upper  integral  is  0  since  the  integrand  is  odd,  and  thus 


m  “  |  (N  »  1) 

We  can  evaluate  the  second  moment  by  writing 


Using  Sterling's  approximation  this  becomes 


(N/ 2) 2  + 


N/2 


i  1 

m  am 

,1  -  m2/ (N/ 2) 2 

N/2-1/2  (1  "  4m/,N) 

N/2 

r  . 

dm  (1  -  4m/N)m 

i 

1  -  m2/(N/2)2^ 

(N-l)/2 

r(N/2  +  m) 
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By  making  approximations 


(l  -  Wto”  -  e-4m  /N 


and 


/  2\^-^  2 

fl  -  4=y  .  e‘2m  /N 

and  extending  the  limits  of  integration  to  infinity,  we  have 


CO 

/ 


2  -2m  /N 

m  e 


m2  =  (N/2) 2  +  0 


/ 


-2m2 /N 


(N/2)2  +  | 


The  variance  is 


m2  -  (m)2  =  (N/2)2  +  N/4  -  (N/2)2  =  N/4 

and  the  RMS  shift  in  the  median  of  any  measurement  from  its  long-term 
average  is  /l/4N  =  . 
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